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■^j- , Abstract 
CO 

fT) ' Field-correlator method is used to calculate nonperturbative dy- 

. namics of quarks in a baryon. General expression for the 3q Green's 

\ function is obtained using Fock-Feynman-Schwinger (world-line) path 

integral formalism, where all dynamics is contained in the 3q Wilson 
loop with spin-field insertions. Using lowest cumulant contribution 
Qh| for Wilson loop one obtains an Y-shaped string with a deep hole at 

Q_|' the string junction position. Using einbein formalism for the quark 

kinetic terms one automatically obtains constituent quark masses, cal- 
culable through the string tension. Resulting effective action for 3q 
plus Y-shaped strings is quantized in the path-integral formalism to 
^ ■ produce two versions of Hamiltonian,one -in the cm. and another in 

the light-cone system. Hyperspherical formalism is used to calculate 
masses and wave functions. Simple estimates in lowest approximation 
yield baryon masses in good agreement with experiment without fit- 
ting parameters. 
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1 Introduction 

Baryons are for a long time an object of an intensive theoretical study [T]-|ll|. 
Both the perturbative dynamics and confinement interaction were considered 
decades ago PQ-[S| and a series of papers of Isgur and collaborators pHE] has 
enlightened the structure of the baryon spectrum in good general agreement 
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with experiment. In those works dynamics was considered as a QCD moti- 
vated and relativistic effects in kinematics have been accounted for. Recently 
a more phenomenological approach based on large A c expansion for baryons 
[12]- [E3 was applied to baryon spectra JTHIIH] and clearly demonstrated the 
most important operators forming the spectrum of 70- plet. 

Summing up the information from the quark-based model one has a pic- 
ture of baryon spectra with basically oscillator-type spectrum, modified by 
presence of spin-dependent forces and other corrections. E.g. the orbital ex- 
citation with AL = 1 "costs" around 0.5 GeV, while the radial one (actually 
two types) amounts to around 0.8 GeV. Moreover, hyperfine splitting, which 
in experiment is large (for A — A system it is 0.3 GeV) is underestimated us- 
ing perturbative forces with a s ~ 0.4 and spin-orbit splitting typically small 
in experiment, needs some special cancellations in theory [T8] . 

Moreover, some states cannot be well explained in the standard quark 
models. A good example is A* (1440), which is too low to be simply a radial 
excitation and moreover its experimental electroproduction amplitudes 19 a 
are in evident conflict with theory |20j . 

This example is probably not unique, and one can notice an interesting 
pattern in "radial" excitations of A, A, A and £: i n a U cases three lowest 
states Mi,M 2 ,M 3 have intervals A, - M 2 - M, - 400 + 500 MeV, A 2 = 
M 3 - Mi ~ 600 - 700 MeV. 

Namely, in 



A(939) 


A(1232) 


A(1116) 


£(1193) 


A(1440) 


A(1600) 


A(1600) 


£(1660) 


A(1710) 


A(1920) 


A(1810) 


£(1770) 








£(1880) 


Ai = 500, 


Ai = 370, 


Ai = 480, 


A x = 470, 


A 2 = 770; 


A 2 = 690; 


A 2 = 700; 


A 2 = 600 



one can see that A 2 corresponds to usual radia 1 excitation, while the energy 
interval Ai cannot be explained in a simple way in standard quark models. 

One can say more about difficulties with the interpretation of the A(1405) 
state, quantitative descriptions of AA transitions etc. [TTH |2*U]. 
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In this situation it sounds reasonable to apply new dynamical approaches, 
which are directly connected to the basic QCD Lagrangian and where all 
approximations can be checked both theoretically and numerically on the 
lattice. 

Here belongs the field correlator method (FCM) started in [2H (for 
a review see (22] and for more dynamical applications (211 )• It is aimed at 
expressing all observables in terms of gauge-invariant field correlators. Its 
use is largely facilitated by recent observation on the lattice [23 I2H] , that the 
lowest bilocal correlator gives dominant contribution to the quark- ant iquark 
forces, while higher correlators contribute around 1%. The use of FCM for 
meson spectra [2ZI |2E| has shown, that gross features of spectra can be calcu- 
lated through only string tension, while fine and hyperfine structure require 
the knowledge of another characteristics of bilocal correlator - the gluon 
correlation length [21] which is known from lattice data [3U] and analytic 
calculations jSH E2] • 

In dynamical applications of the FCM to baryon spectra two different 
schemes are used presently; the relativistic Hamiltonian method (RHM) and 
the method of Dirac orbitals. The first was suggested in |33J EH] and used 
for baryon Regge trajectories in [33 EE] and for magnetic moments in |3*7j . 
The second method was suggested in [38J and exploited to calculate baryon 
magnetic moments in [3T?j . 

Recently an important element was added to the RHM for baryons, 
namely all spin-dependent forces between quarks have been calculated in 
the same Gaussian approximations |40j . To finalize the RHM for baryons 
one still needs to construct the full baryonic Hamiltonian taking into ac- 
count the energy of string motion, nonperturbative self-energy corrections 
|4lJ etc. The present paper is aimed at the fulfilling this task. It contains 
the detailed derivation of the full baryonic Hamiltonian both in the cm. 
and in the light-cone system of coordinates, simple estimates of spectra for 
spin-averaged masses and a preparatory discussion of future explicit detailed 
calculations. 

The paper is organized as follows. In section 2 the 3q Green's function is 
written down using the Fock-Feynman-Schwinger representation [12], in sec- 
tion 3 the cm. and relative coordinates are introduced and the einbein field 
n(t) is introduced which will give rise to the quark constituent mass. In sec- 
tion 4 the resulting effective action is quantized and the full cm. Hamiltonian 
is explicitly written down. 

Section 5 is devoted to the derivation of the light-cone Hamiltonian, its 
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physical interpretation and correspondence to the partonic model. 

In section 6 the construction of the baryon wave function is discussed, 
expansion of its coordinate part into a sum of hyperspherical harmonics and 
analytic estimates of spin-averaged spectra. 

Spin-dependent forces are discussed in section 7 and section 8 is devoted 
to conclusions. 



2 Baryon Green's function 

One can define initial and final states of baryon as a superposition of 3q 

states 

^in,out(%^ \ X^ \ X^ ^) T i n ^ ou fC aia2 ag^-'^P^ i ^X^ \ X^ ^) (1) 

where a, are color indices, while 7, contain both flavour and Dirac indices, 
and a sum over appropriate combinations of these last indices is assumed 
with T as coefficients. Moreover x^\x^ are coordinates of quarks and of 
the string junction respectively. 

The 3q Green's function can be written as 

G 3q (x®\yM) = (tr Y T out l[S^\y^)T in ) (2) 

where we have neglected the quark determinant and defined 

tfY = 7: /_, ^axa^a^bxb^bz (3) 
e —f 

a,i,bi 

and S^(x^\ is the quark Green's function in the external gluonic fields 
(vacuum and perturbative gluon exchanges). For the latter one can use the 
exact Fock-Feynman-Schwinger (FFS) form |2*2*1 l2"4l 14*2*] 

poo rs 

S(x,y) = (m - D) ds(Dz) xy e~ K W z (x,y) exp g J a^F^zij^dr (4) 

where W z is the phase factor along the contour C z (x,y) starting at y and 
finishing at x, which goes along trajectory which is integrated in (Dz) xy , K = 
m 2 s + \ Jo z 2 dr 

rx 

W z (x,y) = Pexpig A^dz^ (5) 
Jy 
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and P is the ordering operator, while 

<v^=(£e,' *b) (6) 

Here E, B are color-electric and color-magnetic fields respectively. The 
average over gluon fields, implied in (J2J) by angular brackets, is convenient 
to perform after the nonabelian Stokes theorem is applied to the product of 
W z . 

Consider to this end the gauge-invariant quantity 

W 3 (x,y) = tr Y flW Zi (x,y) (7) 

i=l 

We can now write the nonabelian Stokes theorem expressing in W Zi 
(e.g. using a general contour gauge) through F^ u (u, x^) = (j>(x^°\u)F liV (u)(j)(u, x^) 
where (f>(x, y) is a parallel transporter, and it is convenient to choose x^ at 
the common point x = x^ = x^ 2 ' = x^ 3 \ 

Making final points also coincident, y = y^ 1 ' = y^ = y^ 3 \ one can use 
the identity 

tr Y 4>a x w 0, y)<Pa 2 b 2 (x, y)<pa- A b- A (x,y) = 1 (8) 

and rewrite W 3 as (the simplest way to derive Q and subsequent equations 
is to choose in the gauge invariant expressions (JTj), the contour gauge, 
where <f) ab = 5 ab ) 

oo (j n \n r 

(W 3 (x,y))=tr Y expY, - 1 r L ((F(l)...F(n)))ds(l)...ds(n). (9) 

,where in F(i) and the surface elements ds(i) the Lorentz indices and coordi- 
nates are omitted, and double angular brackets imply cumulants,as defined in 
|23j . Note that integration in Q is over all three lobes Si, made of contours 
C Zt (x, y) and the string junction trajectory Zy(s), with zy(0) = y, z Y (l) = x. 
The actual form of zy{s) is defined by the minimal action principle and not 
necessarily coincides with the trajectory of the center of mass of 3q system. 
An important specification is needed at this point. 

For this case of 3-lobe loop as well as for several separate loops one can use 
the following gauge-invariant averaging formula, where both field correlators 
are transported to one point x and a, b, c are fundamental color indices 

it?! \ rpf \ \ (tr(F(u,x)F(v,x))} 1 

{F(u,x) ab F(v,x) cd ) = -= {d ad 6 bc - —6 ab 6 cd ). (10) 

Nf - 1 N r 
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Now whenever F(u, x), F(v , x) are on the same lobe, then indices b and c 
coincide and one obtains 

it-,/ \ ti/ (tr F(u,x)F(v,x)) , , 

tr Y (F(u,x) ab F(v,x) bd ) = ± 1 ' / 1 ' ;/ (11) 

where (JBJ) was used. For u, v on different lobes, one instead has 

,t-, \ / \ \ (tr(F(u,x)F(v,x))) , , 

try x) afeJ F>, x) cd ) = iv c (iv c _ ij ( 12 ) 

As the last step in this chapter, one can include the quark spin operator 
u^vF^v into the cluster expansion ©, with the help of relation 



s 



(F^u, x) exp[ig \^ F Xr7 (v, x)ds x <r(v)\) 
1 5 



{exp[ig / F X a(v,x)ds X cr(v)}). (13) 
J s 



Exponentiating the operator F^ v one arrives at the shift operator exp \(s iw gs ( u \ ) 
and finally gets (cf. Eq. (jlj) 

([Wsexpgj^a^l f F^(«®(rW))dr«]> = (^ 3 exp(^F)) = 

8=1 J0 

tryexpE^/ «F(l)...F(n)»dp(l)...dp(n)] (14) 



where we have defined dp(n) = J2i=i dp^ l '(n) 

dp®(n) = dsf nVn {u^) + -cr^Jr^in). (15) 

Here index i = 1,2,3 refers to three lobes Si of the total surface and it is 
understood that whenever F(i) under the cumulant sign ((...)) is multiplied 
by dr^\ it is taken at the point z^'(r^), lying on the quark trajectory z^\r) 
which forms the boundary of the lobe {€). 
Inserting (Jlj), ()14j) into (j2J one obtains 

G 3q (x,y) = 
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" roc 

tr L [T out Y[(m t -D^) R / d Si (Dz^) xy e~ K ^W 3 exp(gaF))T m }. (16) 

i=l ^° 

Here tri is the trace over Lorentz indices, and (m, — D^)r is the value 
of operator (m 8 — Z)W) when acting on the path integral, which was found 
in to be 

{m i -D®) R = m i -ip®. (17) 

Here is the operator of momentum of the quark i. Eqs. f|16|) . ()14|) 
give an exact and most general expression for the 3q Green's function, which 
is however intractable if all field correlators are retained there. 

To simplify we shall use the observation from lattice calculations [23 EH] 
that lowest (Gaussian) correlator gives the dominant contribution (more than 
95%) to the static QQ quark potential. Assuming that situation is similar for 
3Q case and also for light baryons, we now keep in (fT4"j) and only lowest 
cumulant ((FF)) and express it in terms of scalar function D, D\ as in j2I] 

g 2 

—(trF^(u,x)F pX (v,x)) = 
1 d 

= (8np8v\ ~ 6^u P )D(u -v)+ 2^g^~^ u ~ v) p d u \+ perm.)Di(u - v). (18) 

Here we have replaced parallel transporters 0(w, x)(p(x, v) by the straight-line 
transporter 0(w, v), since for the generic situation with \u—v\ ~ T g , \u—x\ ~ 
\v — x\ ~ R, R 3> T g , the former and the latter are equal up to the 
terms 0((T g /R) 2 ). Now in view of JQ, (fTSjl one can write in Gaussian 
approximation 



2 3 

(W 3 exp(gaF)) = expf-^^T J (trF(u)F(v))dp^(u)dp 



g_ 

N C (N C - 1) 
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— - 53 f(tr F{u)F{v))dpV{u)d P «\v)) (19) 

where dp^ is defined in ()15j) . Here (trFF) can be expressed in terms of 
D,Di and one has a closed expression for the term {W 3 exp(gaF)), which 
acts dynamical kernel in the path integral (|16p. 

Now for large sizes of Wilson loop W 3 , such that R {i) > T 9 , one can 
discard D\ and retain D in (I18j) . since only the latter ensures area law (and 
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moreover, lattice data [30J show that D\ C 5). Then the diagonal terms in 
the sum of the exponent in (|19|) can be written as (neglecting spin-dependent 
part for the moment) 

(W 3 ) dmg = exp(-a(S 1 + S 2 + S 3 )) (20) 

where Si is the area of the minimal surface between trajectory of quark 
(z) and trajectory of string junction (Y-trajectory), and we have used the 
relation |2H 122] 

1 r*. 



a = - J d 2 xD(x). (21) 

Let us turn now to nondiagonal terms in f)19|) . Since D(x), Di(x) are 
dying exponentially fast for x > T 9 [Hn]-j32], only region of width T g around 
Y trajectory contributes to these terms, which one can write as 

T/ rp _ ST* TfAij 

Vnondiag.-L v nondiag. 
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Nc{ £^Y } E / (trF^u)F pX (v))ds$(u)ds$(v). (22) 

Separating out time components, u = («4, u), v — (t>4, v), one can write for 
the D contribution 

^nondiag. 



K — x E / D(u® - v^)d(uf - W j*>)d (^^^J d{uf - (23) 

where we have introduced for u, v parallel and transverse components on the 
lobe Si with respect to the lobe Sj. 

Since | U W - = sj(uf - v^) 2 + (uj - u^f + (w? - u { [ y ) 2 is grow- 
ing fast with — t>j_| one can estimate (J2Uj) . (J23*Jl as 

K^L, - Vfi. = C7(r« + + r< s >). (24) 



where rW = zW — z^ y ^ , i.e. the difference of quark and string junction 
coordinates at a given moment of time. Being always smaller than for 

large nevertheless V^ di brings about an interesting cancellation for 
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small r^'. Estimating integrals in ()23|) for small r®, < T g , and using for 
D(x) the Gaussian form, D{x) = D(0) exp(— JU), one has 

9 

V con f = Vdiag. + V non diag. = 7T~/=^tC^2 R^) 2 - (25) 

Z V nI 9 i 

It is clear that for a symmetric configuration = = one has 
Vconf = 0. To study further this cancellation let us take into account that if 
the triangle made of quarks has all angles less than 120° (since string junction 
is at Torricelli point)then the string junction is inside this triangle. In this 
case one can write 



= 2((r«) 2 + (r^) 2 + (r®) 2 ) + r«r^ + r^r^ + r<V 3 > (26) 

i.e., V con f vanishes quadratically in differences of quark distances from the 
string junction. Practically this brings about a strong effective cancellation 
in V con f for 3q system with equal masses at approximately equal distances. 
Numerically and analytically this fact was discovered first in jUj for static 
3Q potential. It was argued there that V non diag brings about smaller slope of 
V(3Q) at small to intermediate distances, as was indeed found on the lattice 
[35j . An assumption that a triangular 3q string configuration is responsible 
for the smaller slope, however, cannot explain it, since that configuration is 
impossible to construct in a gauge-invariant way jl^l- Explicit expressions 
for V^ in general case are given in jl^]- Here and in jlHl a missing in 
factor of —1/2 in front of V non diag is restored. 



3 Gaussian representation for the effective ac- 
tion of quarks and string 

Consider now the exponent of the FFS representation for the 3q Green's 
function 1)160. (|T§j) in the simplified case when i) spin interaction is neglected 
and ii) large distances |Rj| 3> T g are taken into account. 

In this situation one can use the form ()2U|) instead of ([190. and writing 
the exponential term in (fTB*)) as 

G 3q (x,y) =tr L [T out Y[(m i -ip i ) f ds i (Dz^ ) ) xy T in }e' A (27) 

J 
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where A plays the role of effective action, 



A = J2(K i + aS t ). (28) 



i=l 



Our purpose is finally to construct the effective Hamiltonian, considering A 
as an effective action for 3 quarks and the composite string with the string 
junction. To achieve this goal one must i) go over from the proper time 
S{ to the real time integration in the 4c? Euclidean space-time (later on to 
be transformed into Minkowskian time ), ii) to transform the Nambu-Goto 
form of the lobe area Si (see below in ()36|) ) into a quadratic form, as it is 
necessarily done in string theory (since otherwise the path integral (J27j) is not 
properly defined). Both operations are the same as in the qq case, considered 
in and we shall follow closely that procedure. 

The resulting Hamiltonian depends on the choice of the hypersurface, 
and for the qq system both the cm. jMj and light-cone [HI HH] cases were 
considered. 

Below in the next chapter the cm. Hamiltonian will be derived, and to 
this end we choose the hyperplane intersecting all 3 quark trajectories and Y 
trajectory at one common time t, to be considered in the interval < t < T, 
so that quark coordinates are = (t, z' 1 '), and string junction coordinate 
isz( y ) = M y )). 

Now one can make a change of variables, introducing the einbein variable 
[SHEnj fi(t) for a given trajectory zW(r"), < r < s, one defines 

±£vp 2^4") 

so that kinetic term Ki becomes 

K < = £ d «S) + & r (i2(t) + 1)] - (30) 

The transition from the integral over dsidz® to the integral over Dfi®(t) is 
known to have a nonsingular Jacobian (see Appendix A of the second paper 
in for more details and explanations) 

ZVW ~ expM / Jfi 2 (t)dt]dsDz 4 (t) (31) 

£ JO 
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where e ~ 1/A, and A is an ultraviolet cut-off parameter. 
Hence the integrals in (J27|) can be rewritten as 



JJ d Si D 4 z^ -> [] DmD 3 z® (32) 

where the integration measure for .D/ij can be specified further to be |3*4"] : 
Dfj,(t) ~ 11^=1 M 3/2(f \ ■ A s a nex ^ s ^ e P one introduces the cm. and relative 
coordinates 

£=-4^x>(*)* co (*) 

• /3 ^ lZ W+/i 2 z( 2 ) , (3) 1 . /i lZ W - /i 2 z( 2 

^ = V2 ( 2 ^ • (33) 

Here //+ = £f =1 fr. 

From our discussion above it is clear that the time t coincides with the 
fourth component of the cm. coordinate, t = R4, and the whole quark 
kinetic term in (1251) is: 



3 

Here /i^ and //£ will be found below, (|42|) . The area-law term in (}2*5j) can 
be written as follows 

5>^ = a£ Trft I 1 dPiyl{wf)\w'Vy-{wfwT) 2 (35) 
i=i ^i 70 -A) 

where wjp(t,j3) is the i-th string position at the time t and coordinate (3 
along the string and the dot and prime signs have the meaning of the time 
and j3 derivatives respectively. In the spirit of our approach one should take 
the world sheets of the strings corresponding to the minimal area of the sum 
of surfaces between quark trajectories an Y- trajectory of the string junction. 
At this point we make a simplifying approximation [331 EH] that strings at 
any moment t can be represented by pieces of straight lines. In this way one 
disregards string excitations (hybrids) and mixing between these excitations 
and ground state baryons. This can be done for ground states since the mass 
gap for string excitations is around 1 GeV |2~4"j . 




+ i /U+ (t)R 2 + + - Me £ 2 



(34) 
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For higher excited states the mixing should be taken into account analo- 
gously to what was done in meson sector [qP| . 
Thus one writes 

vf{m = zf(t)P + z<P(t)(l-0) (36) 

and time derivatives of w^' in can be replaced using (J5fi)l by time deriva- 
tives of zjf\ z^\ Since also the string junction position is expressed through 
the quark coordinates, the string does not possess the dynamical degrees 
of freedom of its own (in this straight-line approximation). To recover the 
latter one can use background perturbation theory and consider the states 
with 3q and additional valence gluon(s). The latter describes gluonic exci- 
tation of baryon and has its own dynamical degree of freedom. Note that 
this way of systematic description of string excitation is different from the 
ad hoc assumption that string is described by Nambu-Goto action with all 
dynamical string degrees of freedom included, which does not follow from the 
QCD Lagrangian. 

Consider now string-junction trajectory. In line with the whole approach 
one requires that at any given moment z^'(t) occupies the position which 
gives the minimal string energy, i.e. z^\t) should coincide with the Torricelli 
point, giving the minimum of the sum of lengths of 3 strings: 

L = £|z«(t)-z^(t)|, -^- = 0. (37) 
i= i dz\ \t) 

Therefore z^ Y \t) is not an independent dynamical degree of freedom and 
iS Y ^ is expressed in terms of z^\ % = 1, 2, 3. 

Now one can introduce (as is usual in string theory [HT]) the auxiliary 
fields (einbein fields |49j ) to replace the untractable square-root terms in 
(J55J) by quadratic expressions. In this way one writes 

Si = ^[(^ W ) 2 + WW - 2r/,rfr«) + fa) 2 (r«) 2 ]. (38) 

Here Vi(t,(3) > and rji(t,(3) are two einbein fields( which are integrated 
out to yield back the form ()35j1 ). and i"W = z^' — z^ Y \ As a result one has 
for the 3q Green's function 

3 

G 3q (x, y) = J DRDCDrj [] D^D^D^triT out { mi - ip^T^e^. (39) 

i=l 
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4 Quantization of the strings and derivation 
of the 3^-string Hamiltonian 

The action (|2*8j) using (|3~3j) . (|3*5jl. (|3*Sjl and can be written as 



.4 



2 Jo 2 Jo 



< + ^? + £ + / dj g 

2/ii 2 2 Jo 2v, t 



1 r 1 



(40) 



where we have defined Ui = \jv^ and r$ = zW — \ z^'(t) = (t,z^ >). 

As the next step we introduce the cm. coordinate R and Jacobi coordi- 
nates £, 77 as follows [SH] 



4 1} = ^ + 



1/2 



JLf + (/L*l + /i 2 ), 



6 



/X//2 



1/2 



f (2) 



^3 



1/2 



1/2 



;(3) 



,^2(^1 + H2), 
1/2 



(41) 



with the inverse expressions 



Rk = — , m = (4 
^+ i=i 



•(2) .(1^ 



/i!/i 2 



/i(/i! + ^2), 



/'3 



1/2 



( / i 1 4 1) + /U2 4 ZJ -(^i + /i2)ir) (42) 



■(2) 



1/2 



(3> 



In (|4ip . (|42|) the mass /i is arbitrary and drops out in final expressions. 
Using one can rewrite the kinetic part of the action as follows 



t A'. - f *E (g + f ) + + ^ + « 2 )]- (43) 
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The string part of the action can be transformed using (|38|) and integrat- 
ing over r\i to the form 

exp(-^a^) =exp{- C dt\j^ f dA[^((r (!) A + z (y) ) 2 - 
i=i Jo 2 i=i J 



a 2 v 2 



-[(r«/3 + zW)r«] 2 p^) + -^]}. (44) 

At this point it is important to note that z^ Y ' is not a dynamical variable, 
since it is defined by the minimum of the action. 

Taking this minimum at a given moment, one arrives at the definition of 
z^ Y \t) as a Torricelli point, which is to be expressed through the positions 
z«(i); 

z^(t) = f(z 1 \t),z^(t),z^(t)) (45) 

where the function / is defined explicitly in |3 6 3 . Therefore z^ Y \t) is also 
expressed through z^'(t), or through R(t) and r^\t). 

Below the simplified procedure will be used, where one identifies z^ 
with the cm. coordinate R, which is true on average for equal mass quarks. 
Explicit formulas for the general case z^ ^ R are given in Appendix. We are 
now in position to get the final coordinates r/, £ or their linear combinations 
r M = Z W _ ji read from ([41)1 . To this end we replace in (pHj) iP' by R and 
integrate over DR in both expressions ()43|) . (}4"4*j) in the same way, as was 
done in the second paper of( [31], Eqs. (j3Tj) - (j4*U|) ). with the result 



A = f T dt \ {e f — + A + ti*? + £ 2 ) + E f dPiWm+ 



+^- L + ^!((^) 2 - i5 p^-)]}- («) 

The last term on the r.h.s. of ()46|) can be rewritten as Vif3f- — ^ — — and 
disappears when the partial angular momentum Zj vanishes. 

In this case simplifies and using in (j4Tj) instead of t], one gets 
in a standard way the Hamiltonian 

i 3 2 2 i a 2 (r^) 2 

H o = £ EI— + ^ + - + / dA(^i(A) + ^ 1 -^ L )]- (47) 
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One can now apply to (J47|) the minimization procedure to define /ij, z/j from 
the conditions 

= ^ = ^ (48) 

which yields 



I* = ypf + m h v % = &\ r i\- ( 49 ) 

Note that in this case (Zj = 0, i — l,2,3)z/j do not depend on $ and play 
the role of potential. Inserting (|49|) in (|"""7|) one obtains the form well known 
from the standard relativistic quark model (RQM) 1J- 9 j 



3 



iW = £(VP? + ™? + ^l rW D- (so) 



i=l 



Note that ()50|) is valid under assumptions that (1) string-junction 
coincides with the cm.; (2) J2i=i Pi = 0; (3) all angular momenta of quarks 
li, i — 1, 2, 3 are zero, so that only radial part of momentum Pj enters in ()50|) . 
However, in RQM the form (J5U|) is used without the condition (3). As one 
will see in what follows, at nonzero U the Hamiltonian H will be modified, 
and for not large < 4 this modification can be taken into account as a 
string correction AH str ing similarly to the meson case in [3*4*1 \27\. 

We consider now the general case of U > 0. To this end we separate for 
each transverse and longitudinal components as follows (omitting index 
i for a moment) 

1 



r 2 = -{(rr) 2 + (rxr) 2 } (51) 
and correspondingly define transverse and longitudinal momenta 



A = (prr = (52) 



2 _ _ 

2 r I v» \y v\2 



2 _ IP X r ) , , f joo2 /^ 2 l r X r 



f#M«) 2 ^. (53) 



One can now derive the Hamiltonian from (""3"""|) in the usual way 

£i 2 2( yU ,, + / 1 (iA^(^)) 

^r^^/o 1 ^^ (54) 
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This is a general form for any values of li, the limit of Zj — > is obtained 
in (|5U|) . Now we shall derive the opposite limit li — > oo. As in the meson case 
one can argue that in this case /ij <C z/, and one can use the quasiclassical 
method and retain in (J54~j) only the last three terms, expanding them around 
the stationary point at = Tq \ where 



Si)\2 



If 



1/2 

(55) 



Inserting (J55)) back into (|54j) one obtains the following quasiclassical en- 
ergy of the i— th string, = Sf =1 Ei 

i f^r i 1/2 + i r ^ (56) 



and from the stationary point of = 0, one has 




Vi (Pi) = \i-{i-fir i/2 . (57) 

V 7T 

Inserting ()57jl into (|5fijl one obtains finally the energy of rotating string 

E 2 i =2ira<J¥ i , tt = k{k + l). (58) 

This result shows that our general baryonic Hamiltonian indeed admits sim- 
ple rotating string limit at large li, as it was in the case of mesonic Hamilto- 
nian. 

The difference from the mesonic case is however, that for the 3q system 
one should be careful in proper exclusion of the cm. motion and in quantizing 
angular momenta U, which should add up to a total angular momentum 
L = J2i=i h- To do this one should go over from r^',li to the independent 
Jacobi coordinates and momenta r], lg, \ v which add up to L as L = lg + 1,,. 

To accomplish this task, one should express r%' = — using (|4"T|) 
and Appendix through 77 and lj through lg, 1^, and insert it into (p>4"|) which 
makes the whole expression rather complicated and not very tractable. In- 
stead we adopt here another strategy and consider the contribution of string 
rotation as a correction, similarly to the case of mesons [HI], where this ap- 
proach has proved to be successful up to L = 4 j2Zj. Therefore we shall 
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represent the Hamiltonian (J54j) as a sum of unperturbed term H plus a 
string correction AH string , which should work with accuracy better than 5% 
up to I w 3 + 4, 

H = H + AH s tring (59) 

where we have defined 

^ = |g + f)+ El ^+WKr 2 ,r3) (60) 

and 

V conf =aj2\^-^\ = aj2r ll (61) 

i=l i=l 

The string correction AH string is computed to be 

3 l 2 a < r" 1 > 

AHstnn9 = ~ S 2<^>(^ + §<ar l »• (62) 

The total Hamiltonian for the bound 3q system in the cm. coordinates, 
taking into account only valence quarks, can now be written as follows. 

Htot = H + AH s tring + AH cou i + AH se if, + AH S pin (63) 

where H is given in (JSOJ), AH string in AH spin is given in |4Uj , AH cou i, is 
easily computed allowing for perturbative gluon exchanges in W3, resulting 
in a standard expression 

AHml= .^ E ^_ (64 ) 

As to AH self. 1 it was found in [H] to originate from the (aFaF) correlator 
referring to the same quark line. It has the form (|H], Eq. ()36|) ) 

AH self = --Y,^. (65) 
where rji = 1 for light quarks. 
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5 The light-cone quantization of the 3q sys- 
tem. Derivation of the light-cone Hamilto- 
nian 

The general expression of the 3q Green's function allows to calculate Hamil- 
tonian, corresponding to any prescribed hypersurface, with the evolution 
parameter T orthogonal to it, according to the equation (in Euclidean space- 
time) 

BC 

^ = ~HG. (66) 

In the previous section the hypersurface was chosen to be z± = const., 
and the corresponding cm. Hamiltonian was written in (JHBJ). The obtained 
Hamiltonian is a 3q equivalent of the qq cm. Hamiltonian found earlier in 

ma mug. 

The light-cone version of the qq Hamiltonian was derived in JIT) and 
solved numerically in [3%] . 

In this section we shall follow the same technic as in |Uj to obtain the 
3q Hamiltonian on the light cone. To this end one should choose the hyper- 
surface to be the plane with fixed values of z+\ where we use the following 
convention 

ab = a^b^ = afci — a b = a±b± + a + 6_ + a_6 + (67) 
and a± = ^So.. 

The same decomposition of quark coordinate 2® will be used as in ([41)1 . 
but for the light cone (I.e.) coordinates (J67j) . Again for simplicity we shall 
identify i?„ and so that 



rf = zf ~R, = zf - (68) 
Some kinematical properties of I.e. coordinate to be used below are 

= rffr + zP = rffr + Rp, (69) 
r M -r ± , r + _U, ^ 
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= w\ + 2u>_; w± = r±P + R. 

Having this in mind one can directly obtain the I.e. action from PU|) (cf 
[TFj for the equivalent derivation of qq I.e. action) 



-4, 



lc 



2/x< 2 



f(( J R ± + r«) 2 + 



+2(i?_+r_(*))) + ^ / dft 
2 jo 

+^((r?A + R ± ) 2 + 2(r«A + it)) + Vi rgrf>- 

-2vM(f±Pi + R±) r ? + r-)] } , (70) 
where f^, (// =_L, — ) can be expressed in terms of £, 77 using parametrization 



3 3 
i=l i=i 



z (l) 


— Rn = 


f ^3 


Vxi + x 2 


z(2) 


— Rn — 


( x 3 


\Xi + x 2 



1/2 



1/2 . 



•?'2 



\ 1/2 



+ x 2 ) 



1/2 



i (3) - R 



x 2 (xi + x 2 )_ 
/Xi + x 2 \ 



-) 



V x 3 

One can rewrite (|7()|) in the same form as in |U| 

1 r T 



.4 



2 Jo 



Gb + {ai-R_|_ + 2a\R~ + 2a 2 ±R± + 2a 2 



(71) 



(72) 



-2c 1± i? ± + ^[-2c 2i r5!V® + a 3i (r! 



i=i 



.(0 , „ fJQv 1 m * 



2ciirl ; + a 4i (r\ 



]} 



where we have defined 



3 -i 3 

fl i = XX/^ + / v i(0)dP) = Vai 
i=i Jo *=i 



(73) 
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f 1 I 
a 2 ± =Y.(jJ*+ / Vi{f3)f3df3)f { 

02- = £0i + t Ui(j3)PdP)f® = E <*2i 
i=l J o i=l 

cu. = E t dW)wf, c u = t Vi{P)r]idp (74) 
i=1 - Jo 



3 

7 ' - 



C2i= [ d(3fi((3)(3r] h a 3i = m + [ Vi(f3)f3 2 d[3; 
Jo Jo 

a 4i = / d(3(v i r l 2 i +—). (75) 

We now require as in jlTj that transverse velocity should be diagonalized, i.e. 
the mixed term a 2 ± to vanish. This gives conditions on coefficients x$, when 
a 2 ± is expressed in terms of two independent velocities:^ and fj±. 
This immediately yields expressions for xf. 

= + Jo MP)PdP 

Now one can integrate over nf=i Drji in the same way, as it was done for 
the qq system in with the result 

1 pT 3 rX 2/ W\2 2 

A' lc = t dz + { ai (Rl + 2i?_) + E[ / dfi-±±L + + a 3 ,(r?) 2 - 

2 Jo - Jo z/j /ij 



Here we have defined 



(*h2 f r W^2 Jo 



^(/?)(/?-(/?)0 2 ]}- (77) 



(/?)<= fviWPdp/ Cvi{P)dp. (78) 



JO 



Our next step is the integration over DR, which is done in the same way 
as in 07] , and choosing the system where transverse total momentum vanish, 
P_i_ = 0, one obtains 

2 Jo fei 1 • y ° ^ 
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4 (r?) 2 (rS°) a (a u - <i/i 0) » 

where we have defined 

(^ (fc) >= [\(P)(P-(P)i) k dP> (8°) 



Integration over £>i?_ with exp(iPj_ J T R_dz + ) (the exponent appearing 
in the standard way when going from Lagrangian to Hamiltonian represen- 
tation) yields important constraint 8{a\ — P+), i.e. 

ai = P+ (81) 

and integration over DR + is trivial, since A\ c does not depend on R + . 

Before doing calculations for the I.e. Hamiltonian, one should go over to 
the Minkowskian action, which is achieved by replacements 

Mi -> Vi -► -iz/ 4 M , 

ai^-iaf, A^-iA M . (82) 

Omitting the superscript M in what follows one obtains for the Minkowskian 
action 

AM) 1 f T , f /.(i)\2 f 1 cr 2 dp (j) 2 
Ac = o / L i + Wi ) - / (rV) 



2 r / > l ■ ™o*v j_ / / 

■/o [ ^ Jo Vi 

-(zf } ) V± '> ' V r " n '- ■ " } . (S3) 

(r?) 2 (r?) 2 (a l4 - (4 0) > J 

From (|83p one can define in a direct way the I.e. Hamiltonian, writing 
4 M) = / dt + LM, H = £p ± q ± " (84) 



One cannot choose strictly speaking as a set of canonical momenta for 
coordinates since the latter are not independent variables, subject ac- 
cording to (jZ2j) to a condition 

X>r?=0. (85) 

4=1 
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Instead the pair of coordinates £±,r] ± is independent, and one can define 
canonical momenta p±\p±\ as 

(0 1 ^ (V) 1 ^ foa\ 

P± = - P± = ~-r — • 86 

We can use nevertheless = \ gr f^ . Then 

P5 ) = EP^, pS^EpS **,, (87) 

i=i i=i 

where c^, are listed from (J72J) 

Cl?? \xi(xi+x 2 )J ' ° 2,? \x 2 (xi+x 2 )/ ' C3,? 
We are now in the position to use (J84|) and calculate the I.e. Hamiltonian, 

3 m? 1 f±a 2 d(3 (p?) 2 - (p?r«) 2 /(r^ 2 



^ 2/ij 2 Jo ^ 2a 3i 
+ 2(^)2^ + 2(r«)2 G3l a 2 J (2^-a 2l )2 

where defined in (171). (170). 

Momenta p^ are not linearly independent, and from ()85|) expressing 
through p^ one obtains the connection 



fL(p(0 _ Cr« + r% P r 7^ J ) = (90) 
i i "■'>' 

where we have defined 



C = - , (91) 
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n _ (^ 2) >(aii-(^ 0) » + (^ 0) >ait()9)? (Q9 , 

~ ( W\2( I (0)\\ 1 ' 

To understand better the structure of the Hamiltonian (J89|) . consider first 
the limit of heavy quarks m, ^> -y/a, in which case as was shown in j37j , the 
inequality holds fii ^> z/j, % — 1, 2, 3. One has from (fTij) . (J75|) an = a>n = a 3i = 
Hi and the Hamiltonian assumes the form 

f Jmf 1 .ia 2 ^ (P?) 2 - (P^r?) 2 /^) 2 
-n^Q = i o h o / 1 o h 

+ 4^/^^ + MVf±(fWl. (93) 

2(r?) 2 A ^ 2{rf)''i M J 

Introducing the dimensionless quantity yi = (which will be shown to be 
independent of /3 and small, <C 1, one has from (f7H)l and (jHH) 



aj '=iry(^; + ^), £(|r + w) = l (94) 
where F = \ Y,Vi- 

This enables one to expand the mass term in (J93)) around stationary 
points in Xi, and keeping the first order term in y± one obtains 

gS^( M2 < i+2y ) +2M Sfe-i) 2 ^-} w 

where M = Y%=i m i- 

We now define as in [37| the z- component of momenta 

M(x t - ^) = pW (96) 



Keeping now in expansion of (f9"3*J) only leading terms one obtains Hhq 



Vl - with total mass operator 



M 2 = M 2 + 2MY ^—L + M 2 Y ^L- (rf) 2 + 



1=1 



2m 4 fe(rf) 2 
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:P+r W ) 2 a 2 ((rf)^ ' 



M 2 M 2 

One can now define the stationary point of (|9*7j) with respect to yi 



(97) 



(Q) = a ( r l f (r (i) ) 2 = (r (p ) 2 + (i) (Q = P+r® (Q8) 

Inserting y^ = yf*' back into (|97|) one arrives at the familiar nonrelativistic 
expansion 



M 2 = M 2 + 2M^T 



i=l 



'^l + „■<■> 

2m,- 



(99) 



Connection between pW also simplifies for 2/i <S 1, so that ()90|) goes over into 
a simple relation Yh=i — 0, as expected. 

We now turn to the case of light quarks, where relations (|T6j) and (}8T| 
hold with yi nonzero, and observe that three strings contribute to the total 
momentum P + an amount, P+ T = J2i So u i(f3)dl3 which can be significant 
and comparable to that of valence quarks, Yli=i f^i- The numerical value of 
(y) ~ 0.2 obtained in |l8j for a light meson, suggests that a larger value can 
be obtained for P+ r /P+, which can be comparable to the 55% of the energy- 
momentum sum rule, observed in DIS experiment on nucleons. We suggest at 
this point following ^H] that this 55contribution, is mostly due to the string 
contribution P+ r from all Fock components of nucleon, most importantly 
from the ground state strings, and from hybrid baryon excitation, where the 

pstr 

ratio -p— is even larger. 

This point will be elaborated elsewhere [5^] . 

Another topic connected to the I.e. Hamiltonian (|89J) is the whole range 
of dynamical calculations similar to those done for mesons in [3H]- One can 
solve for the eigenfunctions of (|89)l and calculate form-factor and valence 
part of structure function of the baryon-for ground and excited states. To 
compute the full structure function however one needs higher Fock compo- 
nents and first of all lowest hybrid excitations. Here comes in the important 
problem of small x - Regge-type behaviour and connection of t-channel Regge 
poles (including Pomeron) with s-channel summation over baryon resonances 
(primarily hybrid excitation) which is planned to discuss in [5^] . 
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6 Discussion of the cm. wave-function prop- 
erties 



In this chapter we consider possible strategies and first estimates in the de- 
termination of eigenvalues and eigenfunctions of the cm. Hamiltonian (J63j) . 
The latter is the sum of spin-independent part (the first four terms on the 
r.h.s. of (jHSJ)) and AH sp i n , which is calculated in j3U] and has a full relativistic 
Dirac 4 <g> 4 <g> 4 structure. 

At this point one can apply two different approaches in treating AH spin . 
In the most part of this section we shall consider AH spin as a correction which 
should be taken into account in the first order of perturbation theory. This 
is especially consistent for the perturbative part of AH spin , which is known 
for light quarks only to the order 0(a s ). In the next section we also consider 
another strategy, when AH spin , and especially its hyperfine part, is treated 
in a full matrix form. 

We start with the first approach and concentrate on the first term H 
in (|63|). which is given in (|6U|). This term was considered before in and 
numerical solution of the A-type states is presented there including the study 
of Regge trajectories. 

Since H does not depend on spin and isospin and color degrees of freedom 
are already integrated out, one should look for fully symmetric wave functions 
depending on spin variables a, isospin variables r, and coordinates rj. 
Relativistic effects are taken into account in kinematics, where einbein fields 
fii are introduced in (pjUjl .and the latter upon stationary point optimization 
in /ij, yield relativistic energies as in ([50)1 . However it is more advantageous 
to solve (|6U|) . which has nonrelativistic form without square roots, and do 
optimization in /ij for the resulting total mass M (/ij). The accuracy of this 
procedure for mesons was checked in 53, to be around or better than 5%. 
This type of procedure also simplifies calculation of all 4 corrections in ()63|). 
which contain /ij explicitly. 

Hence one can follow the construction of the fully symmetric wave func- 
tion as was done in (HUM], which we slightly simplify and adopt for notations 
used before. Namely Jacobi coordinates 77 (}4~Tj) are chosen to be symmet- 
ric (s) and antisymmetric (a) with respect to interchange indices 1 and 2, 
and belong to the two-dimensional mixed representation of the permutation 
group S3, denoted ip" and ip' while one- dimensional ones are ip s and ip a . 
The same holds true for isospin wave functions 77" , rj' , rf , rf and spin-isospin 
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wave functions ip", <p',p s , p a and finally the full coordinate-spin-isospin wave 
function which should be symmetric in interchange of all 3 indices is 



^(z w , a, r) = + + W + W- 



(100) 



An additional requirement is that (p® and ip^\ i = ft, f, s, a, must belong to 
given total angular momentum L,rriL and total spin S, ms and isospin /, I3. 

Inclusion of AH sp i n helps to construct the wave function as the eigen- 
function of total momentum J,mj. 

Since the construction of spin-isospin functions for 3 quarks is well-known 
PQ-[H], we consicer here only the coordinate part ip^ l \$,T))- As in [2]- [3], 
[3 Eni we shall use the hyperspherical formalism [S3] which has proved to 
be very accurate for the 3q ease, namely the lowest hyperspherical function 
yields eigenvalues with 1 

One can introduce hyperradius p in the following way (note the difference 
from definition in where the case of equal masses was considered). 



P 



E 



i=i p 



(101) 



The coordinate wave function ip(£, 77) can be expanded in an infinite series 
of hyperspherical functions u K (Q) depending on angular variables Q, with 
K- grand angular momentum, K = L, L + 2, L + 4, and v - the set of all 
other quantum numbers, see (S3] for a review, 



P Km 



Writing (jHUjl as 



Pi 



L 2fn 2 



one can reduce the equation hoip 



+ 



dp 2 p dp 



where it was defined 

U u K Up) 



+ 



Eip to a system of equations 
(K + 2) 2 ' 



2pE - 



P' 



Vk 



(102) 



(103) 



(104) 



2P E UkMVk>(p), 



K'v 1 



U 



(105) 
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The confining potential V con f was considered in assuming linear F-type 
form. The analysis of matrix elements (jl(J5j) is also given in [§J and 
here we shall use only the simplest form, namely the so-called hypercentral 
component, which for the equal mass case (/ij = p) is 

U™(p) = 1.118V2ap = 1.58(7p. (106) 

The lowest order equation (|1()5|) for K = K' = was solved numerically in 
[3 121- Below we shall demonstrate a simpler approach which allows to obtain 
eigenvalues of this equation analytically with accuracy of 1% for lowest states. 
To this end we take in ()105|) K = K' (neglecting nondiagonal coupling) and 

making a replacement ip v K (p) = ~~^> reduce equation to the form 

l <fr K {p) 



2/i dp 2 

with 



+ W KK {p)VM = EkuVM (107) 



W KK (p) = bp + ^, 6=.g, d= (K + l)(K + l). (108) 



The eigenvalue Exn can be found using oscillator-well approximation near 
the minimum of Wkk(p) 



dW KK {p) 
dp 

which yields 



p=po = 0, Po = f (109) 



1\ a 2/3 

^Kn = W KK (p )+to[n + -) = -—c n (110) 



2J ^1/3 



where for K = 



3/6^\ 1/3 _ VM 
2 I /i J ' pp\ 



W 00 (p Q ) = - ( — ] ,00=^. (Ill) 



The spectrum cun in (jllUj) corresponds to the so-called "breathing modes", 
when baryon is excited in its p -dependent mode only. 

Finally adding other terms in (|fiOJ) one has for Mx n - the eigenvalue of 

Ho, 

M Kn = -p + E Kn (p) (112) 
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At this stage one defines \x from the stationary point of (|112|) . dK ^ n \/i =l i = , 
which yields 

/2 \ 3/4 [2 \ 3/4 

^° = ( g Cjl ) M KnM = (g Cn ) ■ ( 113 ) 

The total spin-averaged mass of baryon corresponding to the Hamiltonian 
(P)l is 

Mf n = M Kn (fi ) + (AH string ) + (AH coul ) + (AH self ). (114) 

For the lowest states with L — 0,1 one can neglect (AH string ), while the 
other two terms are 



6a _ Xb 1 / 3 

TT/Xq (2/i ) 2 / 3 p 



(AH self ) = -—, (AH eoul ) = - 7K7 ^ (115) 



where (AH se if) is given in jH], while (AH) cou i is in EI] • Here notation 
is used 

3 \ 7T^(7 / 

Now Mxnif^o) is defined in (I113J) and one should choose the only input pa- 
rameters (for light quarks we put all m 8 = 0) a and a s . The string tension a 
is renormalized due to the presence of nondiagonal terms (J23)) and therefore 
is smaller than in the mesonic case (see jlH] for comparison with lattice data 
and more discussion). For simple estimate below we choose a — 0.15 GeV 2 
(the same value as in [Hj) and take a s = 0.4, which near its saturated value 



Results of calculations made according to Eqs. (|112|) - (jll5j) are given in 
Table 1. 



Table 1 

Baryon masses (in GeV) averaged over hyperfine spin splitting for 









a - 


= 0.15 GeV 2 , a s = 


0.4, rrii = 


0. 




State 


M Kn + (AH) sdf 


(AH) cou i 


M tot 

1V1 Kn 


M'°*(exp) 


K = 0, 


n 


= 




1.36 


-0.274 


1.08 


1.08 


K = 0, 


n 


= 1 




2.19 


-0.274 


1.91 


? 


K = 0, 


n 


= 2 




2.9 


-0.274 


2.62 


? 


K = L 




1,71 


= 


1.85 


-0.217 


1.63 


1.6 


K = 2, 


n 


= 




2.23 


-0.186 


2.04 


? 
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As it seen from the Table the calculated spin-averaged mass |(Mjy + 
Ma) = Mqq* agrees well with the experimental average, the same is also true 
for lowest negative party states with K — L + 1, which should be compared 
with states of N and A respectively. 

We also notice that breathing modes (n > 0) have excitation energy 
around 0.8 GeV while orbital excitations K — L — 1 have energy interval 
around 0.5 GeV. 

7 Problem of spin-dependent forces 

We are now turning to the spin-dependent interaction. For the 3q case 
the corresponding nonperturbative and perturbative terms are given in |4Ti] . 
They have been derived under the only assumption of Gaussian dominance, 
i.e. only contribution of the bilocal correlator (represented by scalar func- 
tions D and Di) was retained in @, Gaussian dominance being supported 
by recent lattice data 12S] • The resulting spin-dependent forces have in 
general the form of a product of two 4x4 matrices, one for each interacting 
quark, and this is the most general relativistic spin interaction. 

The expansion in powers of inverse quark mass was not used in and 
for light quarks the spin-dependent interaction is proportional to the terms 
■j^jj- and higher inverse mass terms, where fa are are the same as in (J103|) 
and ()113|) and have the meaning of constituent quark masses, which grow 
with excitation. For the lowest states /i = 0.37 GeV and grow fast with 
increasing K, L and n. 

Now one could use two types of strategy to implement spin-dependent 
forces. 

1). Since all terms in except the last one AH sp i n are diagonal in Dirac 
indices, then one can calculate spin-independent wave-functions and account 
for spin effects calculating matrix elements (AH spin ) KLn . This procedure is 
actually used by most authors, and one can mention two positive moments 
associated with it. First of all in this procedure one treats spin-dependent 
forces as a perturbation, and it should work at least for high enough excita- 
tion, when fafij in the denominator become large. Secondly, the perturbative 
spin-dependent forces are known for light quarks only to the lowest order in 
a s and therefore it is illegitimate to account those terms in higher than first 
order approximation. 

However doing so one immediately comes to a serious contradiction. 
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Namely the theoretical estimates of perturbative hyperfine splitting for a 
reasonable value of a s ~ 0.4 yield values around a hundred of MeV instead 
of 300 MeV for the A — N case [§]. The phenomenological remedy used 
is to take a s ~ 1 and smearing the hyperfine 5-function take the resulting 
potential to higher orders, which was criticized above. 

To resolve this contradiction it is suggested first of all to take into account 
the nonperturbative part of hyperfine interaction, which was derived in jlUl ■ 
It is known to yield a large part of hyperfine splitting in light mesons fZT\ l2"%] 
and may be large also for baryons. Secondly, it is suggested in addition to 
use another strategy discussed below. 

2). In case when spin forces are important, as was discussed in the hy- 
perfine case, one should take into account that the same type matrix element 
which creates hyperfine splitting, also connects lower and higher components 
of Dirac bispinor. Physically this means excitation of negative energy com- 
ponents of quark wave function in baryon, which is also associated with the 
backward in time propagation of quarks. 

Therefore the strong hyperfine splitting implies also strong negative en- 
ergy component excitation, and the solution of the total Hamiltonian (J63j) 
should be sought for in the form 

where a, fij are Dirac bispinor indices and i, k, I refer to the excitation state 
of a given quark. 

This strategy is equivalent to the full relativist 3-body Bethe-Salpeter 
equation, which was studied in the quasipotential form in [SHj. 

Another possible treatment of the same problem was recently initiated in 
(SHUSH], where 3-fold Dirac equations were derived from the QCD Lagrangian 
for the baryon Green's function. 

8 Conclusions. 

We have derived the 3q Hamiltonian both in the cm. and in the I.e. coordi- 
nate systems. It was demonstrated that the cm. Hamiltonian can be written 
conviently as a sum of a main term Hq and four corrections in (|63p. repre- 
senting rotating string energy, Coulomb energy, nonperturbative selfenergy 
correction and spin interaction respectively. The explicit form of all terms is 
given above, except for the last one, published recently in jJUj- 
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The spin-averaged energy levels have been calculated analytically, using 
hyperspherical formalism yielding accuracy around 1% for energy levels in 
linear confining potential. Results for A — N system are in good agreement 
with experiment. The present paper is meant to be a starting point of a new 
treatment of baryons, where all types of forces are derived explicitly from 
the first principles under the only assumption of the Gaussian dominance. 
The spin-dependent forces derived for the first time in its totality in j^U] 
constitute an essential part of this new approach. 



9 Acknowledgements 

The author has benefitted from a good working creative atmosphere in the 
Theory group of the Jefferson Lab. It is a pleasure to thank all members of 
the Theory group, all staff and especially Franz Gross for a kind hospitality. 
Many useful discussions, remarks, suggestions are appreciated during talks and 
contacts with V.Burkert, C.Carlson, R.Edwards, J. Goity, F. Gross, W.Melnit- 
chouk, I.Musatov, A.Radyushkin, D.Richards, W.Roberts, C.Schat, R.Schia- 
villa and J.Tjon. The partial support of RFBR grants 00-02-17836,00-15- 
96736 and INTAS grants 00-110 and 00-366 is acknowledged. The author was 
supported by DOE contract DE-AC05-84ER40150 under which SURA oper- 
ates the TJNAF. 



Appendix 

We start with the definition of the string junction position z^ y ^ which is 
obtained from the minimum condition of the sum 

f>«|=f>«-z( y )| (A.l) 
i=i i=i 

and yields after differentiating in 

3 (i) 3 

EuI)|=£n W = 0. (A.2) 

8=1 I I i=l 

This implies that three unit vectors are at 120° with respect to each 
other, being in one plane. Therefore one can relate positions of quarks z''^ = 
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z 0) _ z u) to rW as follows 



(z fe'))2 = (r « )2 + (r a))2 + | r «|| r 0-)| 5 j ^ j = i ;2 ,3. (A.3) 
One can finally relate r^to the Jacobi coordinates r/, Eq. (j41|) 

r« = 6^ + Ci r/ + <5 (A.4) 

where coefficients 6j, Cj are given in ()41|) and S = R — is to be found by 
solving (jA.2|) . ()A.3|) . Eqs. ()A.3|) are algebraic and allow to find the lengths 
Vi = |rW|,i = 1,2,3 through the quark positions |z^|, therefore r« will be 
assumed to be found explicitly. To find S we place 3 quarks on the plane x, y 
so that quarks 1 and 2 are on the x axis. Assuming the 3q triangle to have 
all angles less than 120°, one can compute both components 8 X , 5 y in terms 
of Vi. Thus the position of the string junction z( y )(xo, 2/o) is found to be 

_ v^nra ( 2r i + r 2K 

Vo — I = , — ^i — , = (A-^J 

/3rl + (2n + r 2 ) 2 J3r 2 + (2n + r 2 ) 2 



where (xi, 0) is the position of the quark 1. 
Similarly one obtains for S x , 5 y 



r2 % TS - r\ ) (2ri + r 2 ) + 3(r 2 - r 3 )r 2 



3r| + (2ri + r 2 ) 



(A.6) 



5, = r3( ; i + r2) " 2rir2 . (A.7) 
x/3 J3r 2 2 + (2n + r 2 ) 2 



One can see that for the symmetric case t\ = r 2 = r 3 the string junction and 
the cm. positions coincide, S x = 5 y = 0. 

Finally we quote for the convenience of the reader the expression of the 
sum (A.l) through the quark positions only, taken from 



j^ Ti = ^[( z d 2 )) 2 + 3 ( Z ( 3 ) - R) 2 + 2V3\z^ x (z< 3 > - R)|] 1/2 . (A.8) 
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